Factorizing Numbers with the Gauss Sum Technique: NMR Implementations 
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Several physics-based algorithms for factorizing large number were recently published. A notable 
recent one by Schleich et al. uses Gauss sums for distinguishing between factors and non-factors. We 
demonstrate two NMR techniques that evaluate Gauss sums and thus implement their algorithm. 
The first one is based on differential excitation of a single spin magnetization by a cascade of RF 
pulses. The second method is based on spatial averaging and selective refocusing of magnetization 
for Gauss sums corresponding to factors. All factors of 16637 and 52882363 are successfully obtained. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 



II. GAUSS SUMS BY DIFFERENTIAL 
EXCITATION 
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Given an integer, factorizing it or confirming it to be 
a prime, is an important problem in network and secu- 
rity systems [lj. Using classical computers with finite 
resources of memory, the factorization time increases ex- 
ponentially with the size of the number. In 1994, Shor 
discovered a quantum algorithm for factorization that 
takes only polynomial time 0. Shor's algorithm requires 
a data register that is large enough to encode the number 
to be factored along with some ancilla qubits. A practi- 
cal demonstration has been carried out by factorizing the 
integer '15' [H, using nuclear magnetic resonance. How- 
ever, quantum computers capable of implementing Shor's 
algorithm for larger number have not been developed so 
far. 



Another physics-based method for factorizing numbers 
has been proposed recently by the group of Schleich 
0) IS HI- The mathematical basis of the technique re- 
lies on the properties of Gauss sums [3] • For the physical 
implementation, it requires an ensemble of two-level sys- 
tems, which accumulate the individual terms. The sum 
is obtained by measuring the expectation value as an en- 
semble average. As specific systems for the implementa- 
tion of this algorithm, atomic systems driven by resonant 
lasers have been proposed [1, 0, El j and an implementa- 
tion b y n uclear spins 1/2 was demonstrated by Mehring 
et al. Hl|. 



Like in Mehring's paper, we also use nuclear spins driven 
by radio-frequency (RF) pulses, but we use two different 
approaches that demonstrate the flexibility of the Gauss- 
sum technique and require fewer RF pulses. 
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A. Principle 

The factorization scheme of Schleich et al. relies on sums 
of the form 
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which are known as Gauss sums, [14j |. Clearly, the series 
adds to infinity if a is integer, and to zero otherwise. 
Schleich et al. used this property for their factorization 
scheme by evaluating the truncated series 
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where (f> m {j) = (2irm y), N is the number to be fac- 
tored, the integer j is trial factor, and M is a truncation 
number. This truncated series adds up to unity if N/j is 
an integer, i.e. if j is a factor of N. In all other cases, 
the sum is a small number whose value depends on the 
trunction number M: 
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1, if N/j = integer, 
<C 1, otherwise. 



Physical systems that can implement this scheme must 
be described by complex numbers. As a first example, 
consider a pendulum with two degrees of freedom. Its 
excitation can be described by an amplitude and a direc- 
tion angle, which can be represented by a single complex 
number. The individual terms of the Gauss sum are real- 
ized by resonant momentum transfers to the pendulum, 
with the direction specified by the phase of the complex 
number. If N/j is integer, the momentum transfers all 
occur in the same direction and therefore keep increasing 
the amplitude of the pendulum. In all other cases, the 
direction will vary and the individual momenta interfere 
destructively. 
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For the experimental realization, we choose a different 
system, which is easier to realize: an ensemble of spins 
/ = 1/2, which is excited by radio-frequency pulses. The 
Hamiltonian describing the interaction with the RF field 
is 



H 



eff _ 



ui {I x cos 4> m (j) + J y sm0 m (j)}, 



(3) 



where I XtV are the components of the spin angular mo- 
mentum operator I, lj is the strength of the RF field, 
and the phase angles (j) m (j) are equal to the phases of 
the corresponding terms in the Gauss sum. 
If the RF field is applied for a duration r, it rotates the 
spins by an angle 9 — lot around an axis in the xy- 
plane, which is oriented at an angle 4> m (j) from the x- 
direction. In close analogy to the mechanical pendulum, 
the effect of the individual rotations adds coherently and 
reaches a maximum if the rotation axes are the same 
(i.e. if j factors N). In all other cases, the orientation 
of the rotation axis is essentially random and the small 
rotations cancel on average. 

Formally, we describe the effect of a single pulse by the 
propagator 



U m (j) = exp(-iH P „{ f r) 



(4) 



as a measure of the overall rotation. This implies that 
the individual rotation angles are small, S < 1, Un- 
der this condition, the individual terms in equation (5) 
approximately commute and the operator product can 
be evaluated also for non-factors. Like in the mechani- 
cal analog, the individual contributions then do not add 
coherently and therefore the effective transverse magne- 
tization remains close to zero. 



B. Experimental results 

In our experiments, we used the 1 H spin of CHCI3 dis- 
solved in Acctone-Dg as our target system. The experi- 
ment was done on a 500 MHz Bruker Avance 11+ spec- 
trometer. Initially, the sample was in thermal equilib- 
rium. The pulse sequence consisted of (M + 1) [8)<j, m 
pulses with the flip angle 9=1° and a delay of 5/xs af- 
ter each pulse, as shown in Figure [T] Since the reference 
frame is resonant with the Larmor precession of the spin, 
the free evolution operators during the delays become 
unit operators. 

At the end of the pulse sequence, we acquired the free in- 
duction decay signal, which is proportional to the trans- 
verse magnetization. For an accurate determination of its 
value, we calculated its Fourier transform and integrated 
over the resonance line. The resulting data were normal- 
ized w.r.t. the reference spectrum (with <fi m = 0, Vm). 
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FIG. 1: Pulse sequence for the differential excitation method. 

As shown in Figure [TJ the Gauss sum is evaluated by 
applying a sequence of M + 1 pulses, all with identical 
flip angle 9 but variable phase cj> m to the spin system. 
The combined effect of such a sequence of pulses is given 
by the propagator 

U(j) = U M ---U Q 


m=M 

where 9 = lot. If j is a factor, <j) m {i) — 2n7r (with n 
integer) and therefore the pulses add coherently. The 
resulting net propagator is then 



170*) =expH(M+l)0/ x ]. 



(6) 



As long as the total rotation angle 9(M+1) remains small 
compared to ir/2, the signal observed in the NMR spec- 
trometer is proportional to the angle and can be taken 
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Factorization of 52882363 
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FIG. 2: Factorization of 52882363 using the differential ex- 
citation method and the truncation number M = 15. 

As a first example, we chose to factorize the number N = 
52882363. The results are shown in Figure [5] for trial 
factors j between 50 and 120. It can be clearly seen 
that the factors are 67, 79, 97, and 103. The visibility 
of the data, i.e. the separation between the factors and 
nonfactors, depends largely on the truncation number 
M. For this example, we used a relatively small value of 
M = 15. For larger values, the relative separation can 
be increased significantly. 



III. GAUSS SUMS BY SPATIAL AVERAGING 

Using the principles outlined above, many similar exper- 
iments can be conceived that achieve the same objective. 
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Echo 




i.e. the magnetization gets back into phase and an echo 
is observed. This holds only for terms where j is a factor 
of N; for the others, the total effect of the sequence of 
M + 1 pulses is not a phase reversal, but only a small 
rotation, which cannot refocus the magnetization. 
In the experimental implementation, we used this scheme 
to determine the factors of the number 16637. Figure 
[4] shows the experimental results for the trial factors j 
between 120 and 140, when the sequence is truncated at 
M = 12. Again, we can clearly distinguish between the 
factors (127 and 131) from the non-factors. 



FIG. 3: Pulse sequence for the spatial averaging method. The 
initial pulse is an RF pulse that creates transverse magneti- 
zation. The Gaussian-shaped second pulse is a magnetic field 
gradient pulse, and the subsequent rectangular pulses are RF 
pulses with flip-angle 6 and phase <f> m . 



As a second example, we choose the case that is sum- 
marized in Figure [3] Here, we also encode the individual 
terms of the Gauss sum in the phase of RF pulses. In this 
case, we choose the flip angle 9 of the individual pulses 
such that the total flip angle adds up to 



(M + 1)9 = 7T 



(7) 



for factors. (More generally, it may be any angle (2n+l)7r 
with n integer). 

This sequence of pulses is applied to a state where trans- 
verse magnetization has been dephased by a magnetic 
field gradient pulse. Starting from thermal equilibrium, 
where the density operator of the system is proportional 
to I z , the excitation pulse creates /^-magnetization. The 
field-gradient pulse turns this into 



I, 



G z 



I T cos a. 



I y sin a z 
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where a z is the positions-dependent dephasing intro- 
duced by the gradient, 



jG z zT, 



(9) 



with 7 the gyromagnetic ratio, G z the field gradient, z 
the coordinate in the direction of the gradient, and T the 
duration of the gradient pulse. The total signal, which is 
the integral of I x over space, vanishes at this point. For 
factors j, the combined effect of the pulse sequence is a 
7r-rotation, which inverts the accumulated phase, inde- 
pendent of the position z [l^] : 



; cos a z 
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I y sma z . (10) 



When the second field gradient pulse acts on this state, 
it adds another phase a z , which is identical to the first. 
The final state is thus 



I r cos a 2 



I y sin a z 



G z 



(I x cos a z + I y sin a z ) cos a z 
— (I y cos a z — I x sin a z ) sin a z 
4, (11) 
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Factorization of 1 6637 
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FIG. 4: Factorization of 16637 using spatial averaging method 
with the truncation number M = 12. 



IV. CONCLUSION 

Interference between different wavepackets is an impor- 
tant feature of quantum mechanics. Schleich et al. have 
used this property to evaluate Gauss sums and shown 
how this can be used to determine if a given number is 
a factor of another number. Their procedure is related 
to the proposal of Clauser and Dowling [13j who used 
optical interferometry for factorization. 
In this paper, we discuss specific physical systems that 
implement such a summation and give experimental ex- 
amples that determine the factors of a given number. In 
particular, we show two NMR experiments that calcu- 
late Gauss sums and apply them to find the factors of 
numbers with 5 and 8 decimal digits. 
The two techniques differ with respect to the conditions 
on the flip angles as well as with respect to the required 
initial condition. We found that the visibility of the 
resulting scans is quite high, even for small truncation 
numbers M ~ 12 — 15. For larger numbers, it might be 
necessary to use larger values of M to obtain sufficient 
contrast between the factors and non-factors. 
The spin system that we use for the implementation 
can be completely specified in terms of a 2-dimensional 
Hilbert space. Accordingly, the dynamics of the system 
can also be described in classical terms and the algo- 
rithm may not be considered a quantum algorithm. This 
fact is also evidenced by the example that we discussed 
(a classical pendulum). In contrast to Shor's algorithm, 
this algorithm is no more efficient than known classical 
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algorithms. 

V. ACKNOWLEDGEMENTS 

This work was initiated by W. Schleich and M. Mchring, 
who also communicated their papers to us before publi- 



cation. We thank Jingfu Zhang for his help in setting up 
the spectrometer and for useful discussions. We grate- 
fully acknowledge financial support from the Alexander 
von Humboldt Foundation, GK-726 Materials and Quan- 
tum Information Processing and the DFG through grant 
numbers Sul92/19-1 and Sul92/ll-l. 



[1] N. Koblitz, A Course in Number Theory and Cryptogra- 
phy ( Springer- Verlag, 1994), 2nd ed. 

[2] P. Shor, in Proceedings of the 35th Annual Symposium 
on Foundations of Computer Science (IEEE Computer 
Society Press, New York, 1994), p. 124. 

[3] L. M. K. Vandersypen, M. Steffen, G. Breyta, C. S. Yan- 
noni, M. H. Sherwood, and I. L. Chuang, Nature 414, 
833 (2001). 

[4] W.Merkel, I. S. Averbukh, B. Girard, G. G. Paulus, and 
W. P. Schleich, Fortsch. Phys. 54, 856 (2006). 

[5] H. Mack, M. Bienert, F. Haug, F. S. Straub, M. Frey- 
berger, and W. P. Schleich, Phys. Stat. Sol. (b) 233, 403 
(2002). 

[6] H. Mack, M. Bienert, F. Haug, F. S. Straub, M. Frey- 
berger, and W. P. Schliech, in Proceedings of the En- 
rico Fermi Summer School, Course CXLVIII, edited by 
P. Mataloni and F. D. Maritini (Elsevier, Amsterdam, 
2002), pp. 369-383. 



[7] W. Merkel, W. P. Schleich, I. S. Averbukh, and B. Girard 
(Unpublished). 

[8] W.Merkel, O. Crasser, F. Haug, E. Lutz, H. Mack, 
M. Freyberger, and W. P. Schleich, Int. J. Mod. Phys. 
B 20, 1893 (2006). 
[9] W. Merkel, H. Mack, E. Lutz, W. P. Schleich, G. G. 
Paulus, and B. Girard (Unpublised). 

[10] W. Merkel, W. P. Schleich, G. G. Paulus, and I. S. Aver- 
bukh (Unpublished). 

[11] M. Mehring, K. Miiller, I. S. Averbukh, W. Merkel, and 
W. P. Schleich, ArXiv:quant-ph/0609174 (2006). 

[12] E. L. Hahn, Phys. Rev. 80, 580 (1950). 

[13] J. F. Clauser and J. P. Dowling, Phys. Rev. A 53, 4587 
(1996). 

[14] Instead of 2, other exponents could also be used. How- 
ever, as in the papers of Schleich et al. we will restrict 
ourselves to 2,which provides efficient randomization. 



